N87- 29944 


EFFECT OF STRUCTURE ON CURRENT AND POTENTIAL DISTRIBUTIONS IN POROUS ELECTRODE* 

Oscar Lanzi and Uziel Landau 
Case Western Reserve University 
Cleveland, Ohio 44106 


Porous electrodes generally contain constricted macropores and localized 
micropores. We have studied the effects of the nacropore constrictions on the 
resistance of a capillary and have developed an analytical model for predicting the 
current distribution in a constricted macropore which directly includes constriction 
effects and does not require an empirical tortuosity parameter. We have also 
investigated the current and concentration distributions in localized micropores and 
have shown that the microporous area is fully accessible to charge and mass transfer 
processes. From these analyses we conclude that the micropores primarily affect the 
kinetics of the interfacial processes by contributing to the interfacial area, while 
the macropores impose ohmic and mass transport limitations through the volume of the 
porous electrode. 


INTRODUCTION 

Porous electrodes are often used in electrochemical reactors because they 
potentially allow electrochemical reactions, which are heterogeneous, to be carried 
out in three dimensions. For example, porous nickel oxide electrodes have been used 
in nickel -hydrogen and nickel-cadmium batteries (in the latter, the cadmium electrode 
is also porous), because the charge and discharge reactions may be made to take place 
throughout the volume of the electrode, thus allowing larger capacities in a given 
cell volume. 

Euler and Nonnemacher [ref.l] followed by Newman and Tobias [ref. 2] developed a 
model to determine the current and concentration distributions within porous 
electrodes under Tafel kinetics, both with and without mass transfer limitations. 
These authors treated the pore structure as a superposition of continuous phases 
without regard to the structural details within the electrode. This approach has 
been been adapted by many workers in models for specific porous electrodes, such as 
nickel oxide electrodes [refs. 3, 4], gas-fed electrodes [refs. 5, 6], and halogen 
electrodes for zinc-halogen batteries [refs. 7-9]. These constitute only a few more 
recent examples; a more extensive review is provided by Newman and Tiedemann [ref 
10 ]. 


The authors cited above followed the precedent of [refs. 1 and 2] by assuming 
the electrode to be a superposition of continuous phases. They did not attempt to 
relate the model parameters, particularly the tortuositj r and specific area, to the 
structural characteristics of the pores. A numerical study based on a random— network 
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model by Kramer and Tomkiewicz [ref. 11] suggests that such relations can be found, 
since their results were found to agree qualitatively with those of a single-pore 
model. The objective of this study is to determine, using both analytical and 
numerical modeling, the significance to the tortuosity and the interfacial area 
parameters of simple one -dimensional models and to relate them to pore structure. 

We will develop a model for prediction of the effect of pore constrictions on 
conduction and diffusion within a non-cylindrical pore geometry, thereby allowing the 
’’tortuosity" parameter to be interpreted in terms of such constrictions. The 
predictions will be compared to a numerical two-dimensional model applied to specific 
constricted pore structures. We also will model current and concentration 
distributions in a micropore embedded in a macropore to demonstrate that the 
microporous interfacial area can readily be penetrated find thus the electrochemical ly 
active area should include it, even when the micropores are quite small. These 
results have important effects on the design of porous electrodes for optimum 
perf ormance . 


SYMBOLS 

Latin 

a Inter facial area/Total pore volume, cm 2 /cm 3 

a^ Interfacial area/Macropore volume, cm 2 /cm :i 

b Tafel parameter, V 

c(z) Local concentration in micropore, mol/cm 3 

c(z) Local concentration in macropore, mol/cm 3 

C ± Integration constant, V/cm 

D Dif fusivity , cm 2 /s 

F Faraday constant, C/eq 

i(z) Interfacial current density, A/cm 2 

i Average interfacial current density. A/cm 2 

avg J 

I'(z) Total current along pore axis, A 

L Length of pore, cm 

R Gas constant, 8.3144 J/mol*K 

S(z) Local cross-sectional area, cm 2 


306 


T Temperature, K 

z Distance into pore from electrode-electrolyte interface, cm 

Greek 

a Transfer coefficient, dimensionless 

A Dimensionless average inter facial current density 

£ Integration variable in Eq. [2], cm. 

n Overpotential, V 

k Electrolyte conductivity, mho/cm 

A Local averaging distance, cm 

<f> Potential, V 

t Constriction factor, dimensionless 

c 

CONSTRICTED PORE STRUCTURES* 

Fig. 1 depicts a schematic model of the structure generally found in porous 
electrodes. There are two scales of pore size: the macropores, which generally 

extend through the volume of the porous electrode structure, and the micropores, 
which are much smaller than the macropores and tend to be confined to localized 
regions such as individual particles. For example, in nickel oxide electrodes pore 
size distribution data [ref. 12] show clearly the existence of 1-pm radius macropores 
and 0.004 to 0.01-pm micropores. We show later that the microporous structure is 
important in determining the active electrochemical area, at least in flooded 
electrodes where the reactants are distributed throughout the volume of the 
electrolyte phase. This section is devoted to the macroporous structure. 

An important feature of the macropore structure is the presence of constrictions 
in the capillary, as can be seen in Fig. 1. An idealized representation of these 
constrictions, on which our study of the constriction effects is based, is given in 
Fig. 2. Such constrictions lead to addditional mass trjinsfer and ohmic resistance 
inside the porous structure. We now pursue a model to describe this constriction 
effect and use it to predict the overall pore resistance. 

Prediction of Pore Resistance 

Following Sides and Tobias [ref. 13] and Lanzi and Savinell [ref. 14], who 
examined constriction effects in the electrolyte phase in bubble layers on gas 
evolving electrodes, we assume that the current in the pore is directed primarily in 
one dimension along the pore axis, and is uniformly distributed in the pore 
cross-section at any point along the axis. We also assume that the microstructure 
can be represented by a segment of length A which is nuch less than the pore length 
L. The model is shown in Fig. 3. With the assumption of a uniform current 
distribution throughout the cross-section S(z), the potential balance takes a simple 
form that is independent of the, pore geometry: 
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d<p/dz = -I'/S(z)k 


( 1 ) 

To integrate this equation it is convenient to introduce the following: 

<f> - (1/A) fl +X f(C) df (2 ) 

The averaged value <f> is essentially the same as f itself if the averaging length A 
is sufficiently small compared to the pore length L. However, derivatives of f and 
<f> with respect to position may differ because of the position-dependent 
cross-sectional area. 


With Eq. (2), Eq. (1) may be averaged over A to give: 

d<<j>>/dz = -<I'/S>/ic ( 3 ) 

Since the averaging length is assumed to be much less than the pore length L, I’ 
may be taken as constant, so that (3) becomes: 


d<<#»>/dz = (-1 ’/k)<1/S> 

For the prediction of pore resistance Eq. (4) is multiplied by A<S>/<S> to 
obtain: 

<*>(z+A)-<Hz) = [-1 'A/k<S>]<SX1/S> ( 5 ) 

The term in brackets represents the resistance of a cylindrical pore segment 
whose length and volume are identical to those of the real pore segment. The 
additional factor <SX1/S> arises specifically from the constrictions in the pore 
structure and hence is a constriction factor. This constriction factor is always 
greater than 1 in noncyl indr i cal pores and thus contributes to the observed 
"tortuosity" factor. The same constriction effect applr.es also to diffusion, so that 
Eq. (5) implies a similar equation for the effect of constrictions on diffusion: 


c(z+A)-c(z) = [-1 ’ A/nFD<S>]<SXl/S> (6) 

We will discuss further the significance of the constriction factor and how it 
relates to pore resistance after showing that this development can be used to 
determine current distributions inside constricted pore structures. Only secondary 
current distributions are modeled, but the analogy between Eqs. (5) and (6) shows 
that this approach works equally well with mass transfer limitations, provided that 
the mass transfer is by diffusion. 

Current Distributions in a Constricted Pore 

Kessler and Alkire [ref. 15] have modeled the secondary and tertiary current 
distribution in through-holes under Tafel kinetics. For constricted pore structures, 
the development of [ref. 15] must be modified to include the constriction effects. 

The current balance for this case is as follows: 

[I ' (z+A)-I ' (z)]/A = dl'/dz = a^<S>i(z) (7) 
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in which the assumption that I* is constant over the length A (see Eq. (4)) 

has been used. Here the quantity represents the interfacial area pre unit of 

macropore volume because we are concentrating on the macropores, which traverse the 
entire electrode length. Combining this with Eq. (4) gives: 

d 2 <<#>>/dz 2 = [(-l/»c)ajji(z)]<S><l/S> (8) 

in which one can observe the presence of the constriction factor already noted in Eq. 
(4). This must be combined with a kinetic equation. For Tafel kinetics: 

dq/di = -d<#>/di = -d<<f>>/di = ±b/i (9) 

in which the Tafel parameter b is RT/aF, and the equivalence between <f> and <<p> 
embodies the assumption that A is much less than L. The ± sign is positive for 
anodic kinetics, negative for cathodic kinetics. As a boundary condition, we assume 
that I* = 0 at the current collector z = L. This is va:.id if the electrode phase is 
a good conductor. Not all porous electrode materials are good conductors, of course, 
but the conductivity of this phase can be improved by, for example, using a sintered 
plate design. With this condition, Eqs. (8) and (9) may be solved analytically to 
give: 


I i I = (C 1 2 (cb/4a M Lr c ) sec 2 [C 1 (L-z)/2b] (10) 

in which t c is the constriction factor <SX1/S>. The integration constant C* can be 
determined from the average current density: 

{C ± L/2b) tan (CiL^b) = t A ( 11 ) 

where 


A = a M^ i avg^ L2/blC ( 12 ) 

The quantity A represents a dimensionless current density analogous to the parameter 
used by Kessler and Alkire [ref. 15] to describe through-hole plating. It is the 
presence of which differentiates this model from that for cylindrical pores. 

We compared this analytical prediction to the results of a two-dimensional 
finite-difference model [ref. 16] which was applied to the constricted capillary 
shapes shown in Figs. 4 and 5. For both geometries the constriction factor is 1.8. 
We investigated the range 0.3<A<1.7. For a 0.05 cm long macropore in 31% KOH 
electrolyte (*c = 0.64 mho/cm [ref. 17]), with a taken as 0.6 (3), this corresponds to 

a M I iavg I ran &i n S from about 3 to 18 A/cm 3 macropore volume. This range is of 

interest because it is here that a significant nonuniformity develops, as will 
shortly be seen. It also corresponds to discharge rates which are likely to be used 
in a nickel oxide electrode (the results, however, are applicable to other electrodes 
for similar values of A.) For the geometry in Fig. 4, the comparison of the 
two-dimensional numerical and one-dimensional analytical results shows good agreement 
even for a nonuniform current distribution; for the largest value of A the current 
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density ratio i(0)/i(L) is about 3 and is predicted to within 7%. Numerical current 
and potential distributions for some values of A are shown in Figs. 6 and 7. 
Similarly, Table 2 and Figs. 8 and 9 show comparisons and distributions for the 
geometry in Fig. 5. Table 2 reveals greater discrepancies than Table 1, and the 
differences are in the opposite direction. The difference between the geometries in 
Figs. 4 and 5 arises from the fact that in these models X is not much less than L, as 
it would likely be in a real porous electrode; thus there are "end effects" at the 
electrode-electrolyte interface because in Fig. 5 more current must pass through the 
constrictions than in Fig. 4. Since the end effects are much smaller in real porous 
electrodes than in the constricted structures in Figs. 4 and 5, they may be averaged 
out, and this leads to an error of about 7% when i(0)/i(L) is about 3. It can be 
seen that the simple model for pore resistance given by Eq. (4) gives a useful 
relation between the pore resistance and the constriction effects without the need 
for an empirical tortuosity parameter. 

The significance of the constriction factor, then, is that it measures the 
contributions of constrictions within the capillary to the pore resistance. In the 
above comparisons, t c could be equated directly to the tortuosity factor. In 

reality, the situation is more complicated. One such effect is suggested by Eq. (7), 
in which the specific area is introduced on the basis of macropore volume. This is 
because it is the macropores that transport mass and charge throughout the porous 
electrode structure; the micropores, as we will see, serve to increase the 
interfacial area but cannot be expected to contribute to the global transport 
throughout the electrode. Thus, the pore-size distribution data in [ref. 12] suggest 
that half of the total pore volume consists of macropores in a typical porous 
electrode; only this fraction can be expected to distribute reactants through the 
electrode. 

Another effect is that of differing pore orientations at various points in the 
electrode. This can be accounted for by using more sophisticated structural models 
than those applied here. One approach is suggested by the random network model of 
Kramer amd Tomkiewicz [ref. 7]. Another is to assume an isotropic cross-linked 
structure whose orientation can be randomized without altering the predicted 
tortuosity. 


MICROPORE PENETRATION 

In the previous development, attention was focused on the resistance or 
tortuosity of the capillary. This section is concerned with the microporous 
component of Fig. 1. The ability of the current and reactant to penetrate these 
micropores is important because this determines what interfacial area is actually 
effective in bring about the electrochemical reaction Jind, therefore, should be used 
in the specific area parameter employed in porous electrode models. 

We examined the current, potential, and concentration distributions in two 
simple micron-size features shown in Figs. 10 and 11. Each contains a micropore of 
equivalent radius 0.1 pm. Micropores may be much smaller than this, but it will be 
seen that this will not affect the results obtained. 

Electrolyte properties employed for the f nite-dif ference calculations 
correspond to the 31% KOH electrolyte used with nickel oxide electrodes. The 
conductivity is 0.64 mho/cm [ref. 17], the transference number of hydroxide ion is 
0.78 [ref. 17], and these data and activity coefficient calculations by Akerlof and 
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Bender [ref. 18] lead to an electrolyte dif fusivi ty of 2.4x10 cm 2 /s. The 

3 

interfacial current density was taken to be +2.4 A/cm total pore volume (compare to 
the previous section where current densities were based on macropore volume). The 
concentration distributions obtained from the finite difference algorithm of [ref. 

16] are shown in Figs. 12 and 13, where c refers to the average concentration in the 
portion of the pore shown in each structure. It can be seen that there is very 
little concentration variation inside the micropore. The same holds for potential 
and current density, for which variations could not he resolved even with four 
significant figures (the concentration variations occur in the fifth place). The 
small concentration nonuniformities that are observed in the graphs arise not from 
micropore depletion, but from the gradient imposed by the macropore to which the 
micropore is connected in each case. The results apply directly to micropores of 0.1 
pm equivalent radius, whereas they are an order of magnitude smaller in actual nickel 
oxide electrodes. However, if the micropore in Fig. 10 or 11 is made 10 times 
smaller to simulate the pore size to be expected in the actual electrode, the 
voltage, current density, and concentration variations within the interior of the 
micropore (excluding the macropore coupling noted above) would only be increased by a 
factor of 10 and would remain insignificant. It follows that nonuniformities do not 
develop in the microporous structure and that the microporous area must be included 
with the electrochemically active area. For example, observed interfacial areas in a 
nickel oxide electrode [ref. 3] are generally two orders of magnitude greater than 
would be expected on the basis of the macropore structure alone. The physical reason 
for this is that the micropores, being much smaller than the macropores, do not carry 
reactant over the length of the electrode. Rather, the}' distribute it locally, on 
the micron scale, and the diffusion and conduction path in the micropore is two to 
three orders of magnitude smaller than in the macropore. This overrides the 

difference in the radii of the individual pores. By the same token, however, these 
micropores do not contribute to the distribution of the reaction throughout the 
electrode, so that when macropore transport limitations set in, the micropores will 
not increase the electrode capacity or the current tliat can be delivered. In 
addition, the microporous paths are far more numerous: even though the micropores 

are individually about seven orders of magnitude smaller in volume, in aggregate they 
occupy about the same volume as the macropores. 

CONCLUSIONS 

Three important conclusions can be drawn from this study concerning the nickel 
oxide electrode and porous electrodes in general. First, the details of the 
macroporous structure are important in determining the ohmic and mass transfer 
resistance within the macropores. In particular, constrictions within the micropore 
structure contribute to the apparent " tortuosity" of the capillaries and this 
contribution can be predicted with good accuracy given a reasonable pore model. 
Second, the macropores alone determine the ohmic and mass transfer limitations in the 
porous structure. These effects are negligible in micropores because these are too 
localized. Third, the interfacial processes occur throughout the microporous area, 
which is thus included in the effective surface area of the pore. Therefore, the 
micropores affect primarily the kinetics of the electrochemical reaction and other 
interfacial processes. Whenever increasing the micropore volume and area fails to 
improve cell delivery or capacity, it is because the limitations are associated with 
macropore ohmic or transport resistance, not with micropore inaccessibility. 
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Table 1. Comparison of one-dimensional analytical predictions (1-dim.) of 

currrent nonuniformity to two-dimensional finite difference results 
(2-dim.) in pore geometry shown in Fig. 4. 


A 

ifOVifU r 1-dim. 1 

ifOVifU T2- 

1.657 

3.286 

3.002 

1.226 

2.532 

2.370 

0.828 

1.938 

1.858 

0.538 

1.565 

1.532 

0.339 

1.336 

1.320 


Table 2. Comparison of one-dimensional analtyical predictions (1-dim.) of 

current nonuniformity to two-dimensional finite difference results 
(2-dim.) in pore geometry shown in Fig. 5. 


A 

i(0)/i(L) n-dim.l 

ifOVifU r 2-dim. 1 

1.425 

2.867 

3.579 

1.043 

2.247 

2.655 

0.731 

1.808 

2.037 

0.492 

1.510 

1.638 

0.317 

1.313 

1.381 



Figure 1- Structure of particulate porous electrodes showing macropores and 
localized micropores. 
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Figure 2. Idealized representation of constricted pore structure. 
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Figure 3 



Constricted pore geometry showing assumed one-dimensional current flow. 
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Figure 5. Pore geometry for constriction effect analysis. This geometry differs from Figure A 
because it has a constriction at the electrode-electrolyte interface. 
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Figure 6. Numerically determined current distributions in the pore ^ 
geometry shown in Figure 4. Numerical values represent A= a^|i^^|L b/<. 
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Figure 7. Numerically determined potential distributions in the pore geometry 
shown in Figure 4. Numerical values represent A. 
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Figure 8. Numerically determined current distributions in the pore geometry 
shown in Figure 5. Numerical values represent A. 
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Figure 9. Numerically determined potential distributions in the pore geometry 
shown in Figure 5* Numerical values represent A. 








0. 00201 


0. 0015 


,j 0.0010 

t 

u 


0. 0005 


0 . 0000 

0.0 0.1 0.2 0.3 0.4 0.5 

Microns from Mocropore 

Figure 12. Concentration variations in the micropore in Figure 10. The 
average concentration in the macropore-micropore structure is c; the variation 
is seen to be quite small. 
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Figure 13. Concentration variations in the micropore in Figure 11. 
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